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mality  must  be  replaced  with  that  of  efficiency.  A solution 
is  said  to  be  efficient  if  it  satisfies  the  constraints 
and  no  other  solution  satisfying  them  scores  as  well  with 
respect  to  all  objective  functions  and  better  with  respect  to 
at  least  one  of  them.  In  the  presented  algorithm,  the  problem 
variables  are  partitioned  into  two  those  whose  coeffi- 

cients in  the  objective  fvinctions  ar«.  all  of  the  same  sign, 
and  the  remainder.  A tree  search  implicit  enumeration 
algorithm  based  on  this  partition  is  developed  and  computa- 
tional results  are  presented. 
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An  algorithm  for  zero'-one  integer  programming 
problems  «ith  more  than  one  objective  functions  is 
developed,  isplemented  and  tested.  For  a 
multiob jective  problem  the  notion  of  optimality  must 
be  replaced  with  that  of  efficiency.  A solution  is 
said  to  be  efficient  if  (1)  it  satisfies  the 
constraints  and  (2)  no  other  solution  satisfying  them 
scores  as  well  with  respect  to  all  objective  functions 
and  better  with  respect  to  at  least  one  of  them.  In 
the  presented  algorithm,  the  problem  variables  are 
partitioned  into  two  sets;  those  whose  coefficients  in 
the  objective  functions  are  all  of  the  same  sign,  and 
the  remainder.  A tree  search  implicit  enumeration 
algorithm  based  on  this  partition  is  developed  and 
computational  results  are  presented. 
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Before  addressing  the  aultiob jective  function  prcbles/ 
the  tree-search  aethod  for  one  objective  function  probleas 
vill  be  reviewed,  nach  work  has  been  done  in  this  area,  and 
different  algorithas  are  described  in  references 
[ 1,2, 3, 4, 5, 6, 7, 8]. 

The  probiea  generally  is  foraulated  as  follows: 


ain 

Z=  C X ♦ 0 X + ...  ♦ 
11  2 2 

C X 

n n 

(1) 

s.t. 

a X +a  X + . . . +a 
i1  1 i2  2 

X >b  i=1,2,. 
in  n i 

. . ,n 

and 

X takes  the  values 

j 

0 or  1 for  all 

j- 

By 

reassigning  subscripts 

and  applying 

suitable 

transformations  on  the  variables  the  problem  can  always  be 
transformed  to  meet  the  following  additional  reguireaents. 

a)  c >0  for  all  j.  (If  some  c <0  then  we  substitute 


b)  c >c  if  l>k. 
1 k 


For  a zero-one  integer  program  there  are  2 candidate 


solutions.  All  these  solutions  are  ordered  in  a diagram  as 
shown  in  Fig.  1 for  n=4. 


Each  node  in  the  graph  of  Fig.  1 represents  a candidate 
solution.  Inside  each  node  there  are  indices  indicating  the 

solution  vith  z =1  for  these  indices  and  x -0  otherwise.  An 

j j 

index  i is  also  associated  with  every  arc.  This  index 

indicates  the  variable  x =1  for  the  node  where  the  arc 

i 

terainates,  and  x -Q  for  the  node  where  the  arc  starts. 

i 


Values  of  all  other  variables  are  the  same  in  both  nodes 


Define  level  i of  the  graph  to  be  the  set  of  nodes  which 
have  i digits  representing  them.  By  convention  level  0 is 
the  level  which  has  only  the  node  0.  It  is  easily  verified 
that  if  there  are  n variables  the  highest  level  will  be 
level  n. 

n 

Note  that  level  i contains  ( ) nodes  and  that  there  is  a 

i 

syaaetry  in  the  structure  of  the  graph  sc  that  the  level 

n/2,  for  n even,  or  the  level  (n/2±1/2) , for  n odd,  have  the 
kiggest  nunber  of  nodes  and  this  nuober  decreases 
symmetrically  as  we  go  from  the  middle  to  the  lowest  and 
highest  levels. 

If  there  exists  a chain  from  a node  Ni  to  a node  Mj, 
then  Ni  is  said  to  be  predecessor  of  node  Nj  and  Nj  is  said 
to  be  successor  of  node  Ni.  All  solutions  are  partially 
ordered  by  the  predecessor-successor  relationship. 

One  solution  is  said  to  "dominate"  another  if  the 
objective  value  assscciated  with  the  first  is  better  than 
that  associated  with  the  second. 


Since  takes  the  value  0 or  1,  the  value  cf  the 

objective  function  Z is  the  sum  of  these  coefficients  c^  for 
which  x^  is  1;  also  since  c^^O  for  all  j,  the  nodes  in 
higher  numbered  levels  represent  worse  (greater)  values  for 

I 

Z than  their  predecessors  do.  Consequently  if  a solution  is 
feasible  or  if  it  is  dominated  by  another  feasible  solution, 
there  is  no  need  to  test  its  successors  since  they  are 
dominated . 


For  exaiple  consider  the  node  1 in  level  1.  Its  Z value 


and  i£  this  sclution  is  feasible  it  dominates  its 


successors  in  level  2,  namely  nodes  12  with  Z=c  -•‘C  , 13  with 

12 

Z^c  '•■c  and  node  14  with  Z=c  4-c  , and  their  successors  in 

13  14 

level  3 (nodes  123,  124  ind  134)  and  in  level  4 (node  1234). 


Another  bounding  relation  appears  from  the  fact  that  the 
objective  function  is  formulated  in  an  increasing  order  of 
the  values  of  the  coefficients  c . So  for  example  if  c <c 

j * ^ 

solution  2 dominates  solution  3 and  solution  24  dominates 

solution  34  even  though  these  solutions  are  not  related  to 
each  other  with  a predecessor-successor  relationship. 

In  references  [7,8]  the  interested  reader  will  find 
example  problems  and  more  details  for  the  one-objective 
problem,  tree-search  type  algorithms. 


t 

1 

I 

h 
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II.  THE  MCLTIOBJECTIVE  PDMCTION  PROBLEM 


When  there  are  more  than  one  objective  function,  the 
notion  of  optima  must  te  replaced  with  that  cf  efficiency. 

A solution  is  said  to  be  "efficient”  if; 

(1) .  It  satisfies  the  constraints,  and 

(2)  . No  other  solution  which  also  satisfies  all  cf  the 

constraints  scores  at  least  as  well  with  respect  tc  all 
criteria  and  better  with  respect  to  at  least  one  of  them. 

A single  objective  implicit  enumeration  problem  will 
have  a unique  optimum  criteria  value,  but  a multicriterion 
problem  can  he  expected  to  have  more  than  one  set  of 
efficient  criteria.  For  example  consider  a problem  with  two 
"minimizing"  objective  functions,  Z and  W.  There  may 

exist  two  solutions  such  that  Z^<Z^  and  Wi>W2  ; in  other 

words  solution  (1)  is  tetter  for  the  objective  function  Z 
and  worse  for  the  objective  function  W.  This  means  that 

both  solutions  must  be  considered  in  the  choice  of  the  final 
solution.  Reference  [9]  gives  an  approach  to  this  type  of 
problem. 

This  thesis  addresses  the  problem  of  finding  all  cf  the 
efficient  solutions,  using  a tree-search  type  algorithm. 
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A.  DIFFEBEMCES  WITH  THE  ONE-OBJECTIVE  FONCIION  PBOBLEH 


i 


Foraulaticn  of  the  probles  for  the  aultiobjective  case 
is  as  follows: 

Bin  2 =c  X +c  X + ...  +0  X 
i i1  1 i2  2 in  n 

s.t.  a X '•■a  X ■«■  . . . -t-a  x >b 
i1  1 i2  2 in  n i 

and  X talces  the  values  0 or  1 for  all  1. 

j 

Clearly  the  constraints  have  the  same  fornulaticn  as 
before;  and  the  only  difference  froo  the  single  objective 
case  is  that  now  there  is  more  than  one  objective  function. 
Because  of  this  difference  the  problem  can  net  be  foraulated 
in  increasing  order  of  magnitute  of  the  coefficients  c 

ij 

To  illustrate  that  consider  the  following  two  objective 
functions: 

Z=3x  ♦Ux  +5x  +... 

12  3 

H=4x  -3x  ♦Sx  +. . . 

12  3 

It  is  clear  that  reordering  N in  an  increasing  order  of 
the  coefficients  c destroys  the  ordering  in  the  objective 

ij 

function  Z. 

The  second  tactic  the  one-objective  function  algorithas 
use  to  reduce  testing,  is  the  formulation  of  the  prcblea  so 
that  c >0  for  all  j.  Unfortunately  no  transformation  can 


i— 1/2/. ..,p 
i= 1,2,...,m 


i 


L. 


i 


12 


i- 


■alee  c 20  for  all  j and  for  all  i;  in  the  ezaaple  above  if 

ij 

one  sabstitates  x =1-7  in  order  to  aalce  c >0  , then  there 

is  an  opposite  effect  in  the  first  objective  function, 

■aking  c <0  . Only  if  c ^0  for  all  i is  this 

ij 

transforaaticn  possible. 

So  in  the  general  case  one  cannot  have  positive  coefficients 
in  all  objective  functions,  and  algorithms  for  the 
■ultiob ject ive  problem  Bust  address  this  greater  generality. 


B.  DEALING  HITH  NEGATIVE  COEFFICIENTS  IN  SOHE  OF  THE 
OBJECTIVE  FOSCTIONS 


Return  now  to  the  graph  of  Fig.  1,  which  has  been 
constructed  for  the  one-objective  case,  and  study  the 
relationship  between  the  nodes  in  the  Bultiob jective 
problea.  Consider  two  nodes  connected  with  an  arc  as  in 
Fig.  2. 

The  solutions  which  are 

associated  with  these  two 
nodes  are: 

Z =c  +c  and 
12  12 

2 =c  *c  *c  =Z  *c 

123  1 2 3 12  3 

It  is  clear  that  the 
relationship  of  these  two 

solutions  depends  only  on  the 

sign  the  coefficient  c has. 

3 

Figure  2 
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If  c >0  then  the  solution  Z doilnates  the  solution  Z 

3 t2  123 


and  if  c <0,  Z is  dominated  by  the  solution  Z 

3 12  123 

Consider  now  two  objective  functions  Z and  H and  suppose 
that  the  coefficient  c >0  (for  the  function  Z) , while  c <0 

13  23 

(for  the  function  W) . No  dominating  relation  can  be 

established  between  the  two  pairs  of  solutions  (Z  ,0  } and 

12  12 

(Z  } since  Z <Z  and  f)  >ii 

123  123  12  123  12  123 

Of  course  if  c was  non-negative  for  all  i the  solution 

i3 

(12)  would  dominate  the  solution  (123)  and  if  the  first  one 

was  feasible  , there  would  be  no  reason  tc  test  the  second 
one. 

The  above  example  easily  can  be  extended  to  the  general  case 
and  the  following  theorem  can  be  established. 

THEOBfi!  1.  In  a multiobjective  problem,  the  solution 

which  is  associated  with  some  node  "a**  in  a 
level  b,  dominates  the  solution  of  some 
successor  of  node  ''a",  node  "aj",  in  the  next 
level  k+1,  if  and  only  if  the  coefficients  c 

ij 

of  the  objective  functions  which  are  associated 
with  the  index  j,  are  nonnegative  for  all  i. 

The  proof  of  this 
theorem  follows  directly 
from  the  above  discussion 
and  will  be  omitted. 

The  notation  is 

illustrated  in  Fig.  3. 
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The  next  theorea#  a direct  result  of  Tbecrem  1,  provides 
useful  insight  into  the  problea. 

THEOBEd  2.  In  a Bultiob jective  problea,  a solution 

associated  with  a node  "a">"ij...",  such  that, 

there  exist  c <0  for  soae  a,  c <0  for  scae  1, 
■i  Ij 

and  generally  there  exists  soae  negative 
coefficient  associated  with  each  one  cf  the 
digits  which  fora  the  node  "a",  cannct  be 
doainated  froa  any  other  solution  in  the  graph, 
so  it  Bust  be  tested. 

PBOOE.  Ihe  proof  will  be  illustrated  with  the  example 
in  Fig.  4. 

Ihe  node  •*123”  can  be 
foraulated  either  from 
node  ”12”  and  the  digit  3 
associated  with  the  arc 
which  connects  the  two 
nodes,  or  frca  nodes  ”13” 
or  ”23”  and  the  digits  2 
or  1 correspondingly.  If 
the  digits  1,2  and  3 are 
all  associated  with  scae 
negative  coefficient,  it 
follows  directly  froa 
theoreB  1 that  no  predecessor  solution  doainates  the 
solution  ”123”,  so  it  aust  be  tested. 

The  problea  will  now  be  reforaulated  and  a tree^search 
algorithB  to  solve  it  will  be  developed. 


C.  FOBHOLAIION  OF  THE  PSOBLEH 


’ The  following  two  transfornations  are  reguired  in  order 

to  foraulate  the  problem  in  a desired  form. 

a)  If  for  a given  j c <0  for  all  i,  then  substitute 

ij 

X =1-y  to  make  the  coefficients  nonnegative. 

j j 

b)  If  c <0  for  some  i and  c >0  for  some  1 and  c >0 

ij  Ij  iJc 

for  all  i,  fcrnulate  the  problem  so  that  j<k.  In  other 

words  shift  the  negative  coefficients  tc  the  left,  by 
renaming  the  variables  or  reassigning  subscripts. 

The  following  notation  will  be  followed  in  the  remainder 
I of  this  paper  when  dealing  with  a problem  with  n variables, 

p objective  functions  and  m constraints. 

SN=  { 1,2,  . . . ,f}  is  the  set  of  the  first  f digits  (f<n) 
for  which  there  exist  at  least  one  negative  coefficient  in 
- some  (but  not  in  all)  objective  functions  associated  with 

them. 

SP={f-«-1  ,f-*'2,. . . ,n}  is  the  set  of  digits  which  are 
associated  with  no  negative  coefficients. 

1 

i 

From  Theorem  2,  it  is  necessary  to  test  all  nodes  which 
have  digits  only  from  the  set  SN,  so  when  representing  the 
set  of  the  solutions  by  a graph,  as  was  done  in  the  one 
oojective  case,  it  is  reasonable  to  consider  these  nodes  as 
a separate  graph.  This  graph  "A"  will  contain  all  nodes 
which  are  combinations  of  the  first  f digits,  so  it  will 
f 

have  2 nodes. 


A graph  **6"  with  nodes  which  are  coabinaticns  of  the  digits 
froB  the  set  SP  can  also  be  constructed. 


IHEOBEtl  3. 


The  set  of  all  possible  solutions  of  an 


integer  zero-one  nultiob jective  function 


problem  can  be  represented  as  the  cartesian 
product  of  the  nodes  of  two  graphs  A and  B. 
Graph  A is  constructed  of  all  coabinaticns  of 
the  digits  from  the  set  SN  and  no  doainating 
relation  exists  between  its  ncdes.  Graph  B is 
constructed  from  all  coabinations  of  the  digits 
froa  the  set  SP,  and  its  nodes  have  the  sane 


predessessor-sucessor  relationship  as  the  nodes 
of  the  graph  which  represents  the  solutions  of 


the  problem  (1 ) . 


f 

Proof.  It  is  clear  that  graph  A has  2 nodes  and  that 
n-f 

graph  B has  2 nodes. 

n 

Their  cartesian  product  is  2 nodes  as  is  expected  fcr  a 

problem  cf  n variables.  No  repetition  of  a node  can  appears 
in  these  products,  since  nodes  from  two  sets  have  no  element 

n 

in  common;  so  the  2 nodes  represent  the  set  cf  all  possible 
solutions  fcr  the  n-variables  problem.  The  fact  that  no 
domination  relation  exists  between  the  nodes  of  the  graph  A 
is  a direct  result  of  Theorem  2.  On  the  other  hand  the 
structure  of  the  graph  B is  analogous  of  the  structure  of 
the  graphs  fcr  the  one-objective  function  problems,  because 
no  negative  coefficients  are  associated  with  the  nodes  of 
this  graph. 

The  cartesian  product  set,  which  is  the  set  of  all 

f 

nodes,  can  be  partiticnad  into  2 subsets,  each  the  product 


of  a node  of  graph  k with  all  of  the  nodes  of  graph  B.  The 
snccession  graph  over  the  resulting  set  contains  only  arcs 
associated  with  all  nonnegative  coefficients;  noraal 
dcainance  tactics  aay  be  employed. 

This  consideration  of  the  problem  by  two  graphs  is  very 
interesting  and  very  helpful  because  it  links  the 
multiob jective  problem  with  the  cne  objective  case. 

This  is  so  because,  as  is  apparent  from  Theorem  2,  in  any 
type  of  algorithm,  all  nodes  of  graph  A must  be  tested;  on 
the  other  hand  as  long  as  graph  £ has  the  structure  cf  the 
graph  fcr  the  one-objective  function  problem,  all  the 
research  which  has  been  done  in  this  area  can  be  used  for 
the  multiob jective  problem,  keeping  in  mind  that  here  the 
coefficients  can  not  be  aranged  in  increasing  order. 

These  two  graphs  provide  an  indication  of  the  sise  of 
the  problem.  Of  course  in  the  worst  case,  it  might  be 

n 

necessary  to  test  all  the  2 solutions,  but  even  in  the  best 
case  it  is  necessary  to  test  all  the  solutions  which  are 

f 

associated  with  the  nodes  of  the  graph  A,  e.g.  2 solutions. 
Normally  f is  a small  number,  because  the  objective 
functions  have  small  inclination  between  each  other. 
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III.  SCiLl  £H0£^£M^  in  £IM£L£^ 


Consider  the  following  exaaple  sulticriterion  prcbles, 
expressed  in  the  canonical  fora  set  out  in  the  preceding 
section: 


subject  to: 


Z1*  X ♦2x  ♦2x  -fix  ♦Ux 

1 2 3 4 5 

Z2»-2x  - X ♦2x  ♦ X ♦ax 

1 2 3 4 5 


x+x+x+x+x>1 

1 2 3 4 5 

-X  ♦ax  ♦2x  ♦zx  -ax  io 

1 2 3 4 5 

X - X ♦2x  - X ♦ X >0 

1 2 3 4 5 


( 3 ) 


X takes  the  value  0 or  1 for  all  j 

j 


A.  TOTAL  SNQHEBATION 


Firsts  all  possible  2 -a2  solutions  will  be  explicitly 

enuaerated  and  soae  statistics  will  be  calculated  in  order 
to  provide  an  indication  of  the  redundancy  obtained  with  the 
tree-search  algoritha  in  the  required  work.  For  this 
purpose  Table  1 has  been  constructed,  where  first  each 
constraint  is  checked  for  feasibility;  if  one  constraint  is 
not  satisfied  there  is  no  reason  to  proceed  further.  If  the 
solution  is  feasible,  then  its  value  is  calculated  and 
stored  in  the  proper  column.  In  the  last  coluan  are  kept 
the  current  efficient  solutions. 


J 


node 

solution 

feasibility  1 

— ^ , — ! 

value 

e£f lent 
sltns 


TABLE  1 


20 


Note  that  in  this  problen  the  solaticn  (12)  was  found 
ver;  early  and  since  this  solution  dominates  all  other 
solutions  except  the  solution  (3) , the  number  of  comparisons 
for  efficiency  was  reduced  to  24.  Of  course  this  is  net  the 
typical  case.  Figure  5 illustrates  the  bounding 
relationships  between  the  solutions,  and  the  reader  easily 
can  verify  that  generally  more  comparisons  are  required  in 
order  to  obtain  the  set  of  efficient  solutions  ( (3)  and 
(12)  } for  the  example  problem. 


E.  THE  TBEE-SEARCH  PECCEDDBE 


Consider  now  a tree-search  type  algorithm  in  order  to 
reduce  the  required  amount  of  work.  Following  the  notation 
introduced  for  theorem  3: 


SN={  1,2  } 


SP={  3,4,5  ) . 


As  theorem  3 states,  the  set  of  all  possible  solutions 
will  be  the  cartesian  product  of  the  two  graphs  A and  B of 
Fig.  6. 


* \ 

/ \ 

» \ 


3 1 fl 


4^0 


X ’X‘! 


.1  ./ 


• / 
• / 


0 R A PH  A 


GRAPH  B 


Figure  6 


As  has  been  seen,  all  solutions  cf  graph  A must  be 
tasted  in  combination  with  the  nodes  of  graph  B. 

Now  a logical  seguence  to  visit  the  nodes  cf  the  graphs  is 
as  follows: 

First  a node  from  the  graph  A is  selected  as  the  rcot  of 
the  graph  which  is  constructed  from  the  combination  of  this 
node  with  the  nodes  of  graph  B.  Then  this  new  graph  is 
searched,  testing  these  nodes  which  are  not  bounded.  For 
example  consider  the  composite  graph  formed  from  graph  B and 
node  (2)  cf  graph  A.  First  solution  (2)  is  tested.  If 
necessary,  its  immediate  successors  (23)  , (24)  and  (25)  are 
tested;  some  of  their  successors  may  also  require  testing, 
and  so  on. 

To  reduce  testing  it  is  necessary  to  determine  which 
nodes  are  dominated.  The  following  rules  apply: 

BOLE  1.  If  a node  is  feasible  there  is  no  need  to  test 
its  successors  in  higher  levels. 

For  example  if  node  (24)  is  feasible  it  is  unnecessary 
to  test  nodes  (234),  (235)  and  (2345). 

ROLE  2.  If  a solution  is  dominated  by  one  of  rhe 
current  efficient  solutions  ( i.e.  solutions 
which  up  to  this  point  have  been  found  to  be 
feasible  and  not  dominated  ) then  there  is  nc  need 
to  test  its  feasibility  or  tc  examine  its 
successors. 

Let  us  apply  now  the  preceding  to  solve  the  example 
problem.  The  procedure  is  illustrated  in  Table  2. 
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•ibility 

va  lue 

current 
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sins 

excluded 

nodes 

a 

1.1 

c 

0 

n 

3 

y 

y 

y 

(2,  2) 

(2,  2) 

34,35,345 

4 

y 

y 

n 

(3,  1) 

5 

(4,  3) 

45 

1 

y 

n 

(1/-2) 

13 

y 

y 

y 

(3,  0) 

(3,  0)  ♦ 

134,135,1345 

14 

y 

y 

y 

(4,-1) 

(4,-1)  * 

145 

15 

(5,  1) 

2 

y 

y 

n 

(2,-1) 

23 

(4,  1) 

234,235,2345 

24 

(5,  0) 

245 

25 

(6,  2) 

12 

y 

y 

y 

(3,-3) 

(3,-3) 

123,124,  125.  1234 
1235,12345,1245 

total 

8 

7 

6 

TABLE  2 


First  to  te  examined  are  the  combinations  of  the  node 
(0)  of  graph  A with  all  nodes  of  graph  B,  testing  for 
feasibility  until  finding  the  first  feasible  solution.  Node 
(0)  is  not  feasible  so  node  (3)  is  tested  (see  also  Fig.  6)  . 
This  node  is  feasible  and  according  the  Rule  1,  it  is 
unnecessary  to  test  its  successor  nodes  (34)  , (35)  and 
(345) . From  this  point  on  there  is  a current  efficient 
solution,  sc  when  node  (4)  is  examined  first  its  value 
(3,  1 } is  calculated  and  it  is  determined  if  it  is  bounded 
from  the  current  efficient  solution.  Because  it  is  not,  its 
feasibility  is  examined.  Since  node  (4)  is  nor  feasible  and 
is  not  dominated,  neither  Rule  1 nor  Rule  2 can  be  applied 
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to  exclude  its  successors  from  testing.  The  next  node  (5), 
which  has  a value  ( h,  3 },  is  dominated  by  one  of  the 
current  efficient  values  { 2,  2 ),  so  it  is  net  necessary  to 
test  its  feasibility  and  its  successor  node  45  can  be 
ignored. 

Proceeding  in  this  manner  Table  2 is  completed,  yielding 
the  following  statistics; 

i.  Number  of  examined  solutions:  13 

ii.  Constraints  examined  for  feasibility:  8+7+6=21 

iii.  Calculated  values  of  the  obj.  functions:  2*12=24 

iv.  Compared  solutions  for  efficiency:  19*2=38 

Compared  with  the  results  from  Table  1,  here  only  4051  as 
many  solutions  were  examined,  only  2H%  of  the  constraints 
for  feasibility  were  calculated,  and  67S  as  many  values  of 
the  objective  functions  were  calculated.  Note  that  here, 
the  solution  (12)  which  dominates  most  of  the  other 
solutions,  was  the  last  one  which  has  been  tested,  while 
before  it  was  tested  very  early;  nevertheless  the  number  of 
comparisons  made  was  significantly  .reduced . 

Here  are  some  more  rules  which  improve  the  efficiency  of 
the  algorithm  in  small  scale  problems. 

ROLE  3.  Consider  a graph  B containing  a feasible 
node  (a)  which  dominates  another  node  (b)  . If 
some  solution  (ma)  is  feasible,  there  is  no  need 
to  test  the  node  (mb) , since  it  is  dominated  from 
the  node  (ma)  . 

In  the  example  node  3 dominates  node  5 . Since 
node  13  was  found  to  be  feasible  it  was  possible  to  avoid 
testing  node  15  . 

Single  objective  implicit  enumeration  (ref  [7])  ccmmcnly 


takes  advantage  of  the  following  two  observations;  since 
they  deal  strictly  with  the  constraints,  they  are  directly 
applicable  tc  aulticriterion  implicit  ennmeration  as  well. 

Observasion  1 . Consider  a node  (a)  with  z =1  if  j is 

j 

in  (a)  and  otherwise,  ill  the  successors  of 

node  (a)  must  have  ^ * these 

variables  are  fixed  for  the  successors  of 

node  (a)  ; all  others  are  said  to  be  free 
variables  as  they  may  take  either  c£  the  values  0 
or  1.  If  (a)  is  not  feasible,  it  is  possible 
that  there  are  not  enougph  free  variables  left  to 
satisfy  a given  constraint. 

For  example  assume  that  a given  constraint  is 
-X  -X  +x  +x  >1  and  node  (12)  is  under  consideration.  In 

12  3 4 

this  case  even  with  x =x  =1  the  constraint  is  still  not 

3 4 

satisfied.  Nhen  this  happens,  there  is  no  need  to  test  the 
successors  of  node  (12) . 

Observasion  2.  When  a subset  of  variables  is  fixed, 
then  a given  constraint  may  force  some  other 
variables  to  be  fixed  also. 

For  example  let  a constraint  be  2x  >x  -x  <0  in  an 

12  3 

n-variable  problem.  Consider  node  1.  In  crder  to  satisfy 

the  given  constraint,  all  the  successors  of  node  1 must  have 
x^  = 1 and  Thus  there  is  no  need  to  test  nodes  (12)  and 

(13),  and  among  their  successors  , there  is  need  to  test 
only  node  (123)  and  its  successors. 


IT.  LARGE  SCALE  PROBiEfi§ 
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As  long  as  ths  problem  does  not  have  too  many  variables, 
one  can  easily  constract  the  graphs  A and  B and  keep  track 
of  the  solutions  which  oust  be  tested  or  not.  But  for  0-1 

of  possible  solutions  grows 

(n) . Ihus 


integer  problems  the  number 

n 

) with  the  number  of  variables 
5 

there  are  2 =32  candidate  soluticns,  but  for  n=10 


exponentially  ( 2 
for  n*5 

there  are  more  than  one  thousand  and  for  n=30  more 
billion. 


than  one 


A.  AN  ACCITIVE  ALGORITHM 


From  the  above  discussion  it  follows  that  for  large 
problems,  it  is  necessary  to  use  a procedure  to  generate 
those  nodes  (or  soluticns)  , and  only  those,  which  must  be 
tested. 

The  structure  cf  the  graphs  A and  B,  and  the  nature  of 
the  problem,  suggests  an  algorithm  of  additive  and/or 
recursive  type.  In  the  example  illustrated  in  the  previous 
section,  the  procedure  followed  was  to  test  a node  and  if 
one  of  the  bounding  rules  held,  to  exclude  from  testing  a 
set  of  successor  nodes. 

Here  the  procedure  is  slightly  changed.  A list  of  nodes  to 
be  checked  is  maintained,  and  after  testing  a given  node,  if 
none  of  the  bounding  rules  apply,  the  successors  of  this 
node  in  the  next  level  only,  are  added  tc  the  list.  The 
nodes  at  highest  levels  are  not  added  since,  if  it  is 
reguired  for  them  to  be  tested,  they  will  be  generated  when 
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their  imoedlate  predecessors  are  investigated,  k convenient 
way  to  Iceep  track  of  the  nodes  which  oust  fce  generated,  in 
order  to  be  protected  from  duplications,  is  as  follows: 

Consider  two  nodes  with  the  property  that  the 
designation  cf  the  second  is  the  designation  of  the  first 
plus  an  additional  index  larger  than  the  larger  index  cf  the 
first.  The  second  node  is  said  to  be  a direct  lexicographic 
successor  of  the  first.  The  "successors'*  of  a node  include 
its  direct  lexicographic  successors,  their  direct 
successors,  and  so  on.  All  solutions  are  partially  ordered 
by  this  relationship.  For  the  graphs  which  are  considered 
as  the  product  of  a node  of  graph  A with  all  nodes  of  graph 
B,  each  node  dominates  its  lexicographic  successors.  Links 
of  the  direct  lexicographic  succession  in  figures  1 and  6 
are  shown  as  solid  lines;  they  constitute  a tree  rooted  at  0 
and  spanning  all  the  nodes  of  the  graph. 


The  above  technique  permits  the  calculation  of  the 
values  of  the  objective  functions  and  the  values  cf  the 
constraints  by  the  following  recursive  equations: 

Zo  (i)  =Zo  (i-1)  ♦C(j) 

( 4 ) 

Zc  <i)=Zc  (i-1)  +A(j) 


Where;  Zo:  denotes  the  vector  of  objective  functions 
i:  denotes  the  level  of  the  graph 
C(j):  denotes  the  vector  cf  coefficients  cf  the 
objective  functions  which  are  associated 
with  X 

j 


2c:  denotes  the  vector  of  constraint  values. 
A(j):  denotes  the  vector  cf  the  coeff icierts  of 
the  constraints  which  are  associated  with 

X 

j 
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As  an  illustration,  in  the  exasple  prcblea  froi  the 
previous  section,  the  above  values  for  the  node  23  in  level 
2 are  as  follows: 

Zo(2)*{  4,  1 } 

Zc(2)»£  2,  5,  1 ) 

Now  to  calculate  these  values  in  the  next  level  3 for 
the  nodes  generated  from  node  23,  nodes  234  and  235,  it  is 
only  necessary  to  add  the  proper  coefficients  in  the  values 
of  the  previous  level.  For  node  234,  C ( j) (4) « (3,  i)  and 
A(j)=A(4)  = (1,  2,  -1)  . 

Thus : 


Zc(3)=  Zo(2) +C(4)  = {(4  + 3)  , (1  + 1)}  = (7,  2)  and 
ZC(3)=  Zc(2)+A(4)  = {(2  + 1)  , (5+2)  , (1-1)}  = (3,  7,  0 ) 

Analogously  for  the  node  235: 

Zo  (3)  = ( 8,  4)  and 
Zc|3)  = (3,  2,  2) 

The  following  notation  is  introduced  tc  help  in  the 
foraulaticn  of  a step  by  step  algorithm  which  employs  these 
techniques. 

SMT1,  SNT2  = The  sets  of  nodes  to  be  tested  in  graphs  A 
and  B respectively. 

SES  = The  set  of  currently  efficient  solutions. 

SOB  = Solution  on  hand. 


Now  the  algorithm  can  be  formulated  as  fellows: 


STEP  0:  (initializations) . 

SES^eapty:  SHT2«eBpty; 

SHT1s(  di  d is  a node  froa  graph  A ) 


STEP  1: 


is  soae  node  in  SMT1 


If  SHTI^eapty  then  stop 
SliT1«SMTl-ak  where  ak  i 
whose  the  last  digit  is  k. 
SliT2*{  ak  } 


STEP  2:  (pick  the  SOH) 

If  SNT2=empty  go  to  Step  1 

SNT2sSHT2-ak  where  ak  is  soae  node  in  SNT2 

STEP  3:  (check  for  doainance) 

i.  SOH=  ak 

ii.  Calculate  Zo  for  SOH  using  ( 4 ) . 

iii.  If  SES  is  empty  go  to  step  4. 

iv.  If  SOH  is  bounded  by  soae  sclution  in  SES 
then  SNT2=SNT2-alk  for  all  l=k-1  ,k-2, . . . ,f and 
gc  to  Step  2. 

STEP  4:  (check  for  feasibility) . 

i.  Calculate  Zc  for  SOH  using  ( 4 ) . 

ii.  If  SOH  is  not  feasible  gc  to  Step  5. 

iii.  Put  Zo  for  SOH  in  SES. 

iw.  Elioinate  froa  SES  all  these  solutions 
which  are  bounded  froa  SOH. 

▼ . SNT2=SNT2-alk  for  all  l=k- 1,k-2, . . . , f+1. 

Go  to  step  2. 

STEP  5:  (Generate  next  level  successors) 

SIIT2-S1IT20  {akjl  j»k+1,k+2,.  . . ,n  if  k>f 

j=f+1,f+2,. ..,n  if  k<f) . 

Go  to  step  2. 


I 


Table  3 suaaarizes  the  application  of  this  algocitha  to 
the  previously  used  exaaple. 


STEP  0:  SES=eBpty;  SNT1»{0,1,2,12) ; SNT2=eBpty 
STEP  1:  SMT1>  {1,2,12}  (Mode  ak  is  node  0). 
SMT2«{  0 }. 

STEP  2;  SHT2  is  not  eapty  so  SMT2=SNT2-0=eBpty . 


STEP  3: 


SOH»0 


ii.  Zo(0)s(0,0) 

iii.  SES  is  eapty  so  go  to  Step  4. 

STEP  4:  i.  Zc(0)  -(0,0,0) 

ii.  Since  6^(1, 0,0),  SOH  is  not  feasible. 
Go  to  Step  5 . 

STEP  5:  He  have  k^O  and  £^2  so  SNT2»  (3,4,5} 

Go  to  Step  2. 


STEP  2:  SNT2=SMT2-(3)  »{4,5}  (He  examine  solution  3) 

STEP  3:  i.  SOH=(3) 

ii.  Zo(1)  =Zo  (0) +C(3)  = (0-*-2,  0+2)*(2,2). 
iii.  SES^enpty  so  go  to  Step  4. 

STEP  4:  i.  Zc(1)=Zc(0)+A(3)=(1,2,2)  . 

ii.  Since  B=(1,0,0),  SOH  is  feasible, 
iii.  SES={  (2,2)  }. 

Go  to  Step  2. 


STEP  2:  SNT2»  { 5 } (ve  examine  now  node  (4)  ). 
STEP  3:  i.  SOH*(4). 

ii.  Zo  (1)=  (0+3,0  + 1)  = ( 3,1  ). 

iv.  SOH  is  not  bounded. 

STEP  4:  i.  Zc  (1)  = (0+1, 0+2, 0-1)  = (1, 2,-1)  . 

ii.  SOH  is  not  feasible.  Go  to  Step  5. 
STEP  5:  He  have  k=4  so  SNT2=SNT2U{  45  } = { 5,45  }. 
Go  to  Step  2. 


STEP  2:  SMT23{  45  } (He  test  no.w  ncde  5). 
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STEP  3 


: i.  S0H»(5) 

ii.  Zo(1)  = ( 4,  3 ) 

iv,  SOH  is  bounded  fron  the  solution  (2,2)  . 

SNT2*S1IT2- (4  5)  =e»pty  ( alc=5  and  l=)c-1  = 4 ). 

Go  to  Step  2. 

STEP  2:  SNT2*e«pty  so  go  to  Step  1. 

Pros  this  point  the  reader  should  not  have  any 
difficulty  following  the  way  Table  3 has  been  filled  in. 

Mote  that  in  the  construction  of  Table  3,  in  order  to 
calculate  the  values  of  Zo  and  Zc  for  a node  with  sore  than 
two  digits,  two  or  sore  values  (if  they  appear  in  the  table) 
are  added.  So  for  exaaple,  to  calculate  the  value  of  Zc  for 
the  node  (23) , the  corresponding  values  of  Zc  for  the  nodes 
(2)  and  (3)  are  added.  The  calculation  of  the  value  of  Zo 
is  analogous. 

So,  for  the  nodes  for  which  the  required  infcraation  already 
appears  in  the  table,  the  equations  (4)  can  be  replaced  by: 

Zc  »Zc  -t-Zc  and 
aj  a j 

(5) 

Zo  =ZC  ♦Zo 
aj  a j 

where  aj  is  the  node  with  last  digit:  j and  the  rest  of 
the  digits  in  the  string  a.  Note  that  j can  also  be 
considered  as  a string  of  digits  and  that  the  digits  which 
fors  the  node  can  be  partitioned  in  aore  than  two 
substrings,  for  which  the  values  of  the  corresponding  nodes 
Bust  be  added  to  calculate  the  values  for  the  exaained  node. 


E.  laPIEtlEIAIIOIl  USING  THE  COHPOTEfi 


[ 

r 

j 

• » 


In  this  chapter  the  structure  a cosputer  progras  eust 
have  to  solve  the  aultic^ jective  problen  will  be  ezaained. 


In  the  previous  sections  the  values  cf  Zc  and  Zc  were 
calculated  recursively  using  the  equations  (4) , and  the  next 
level  successors  of  each  node  were  lezicographicaly 
generated  by  successively  concatenating  tc  the  node  all 
digits  which  are  greater  than  its  last  digit.  Norking  with 
the  coaputer  it  is  important  to  keep  in  storage  only  that 
inforaation  which  ,is  required  to  proceed  with  the  following 
steps.  There  are  two  approaches  to  search  a graph  and  to 
generate  its  nodes. 

One  approach  is  to  search  the  graph  level  ky  level.  First 
all  the  inforaation  which  corresponds  to  the  level  1 is 
stored.  Fron  this  inforaation  for  the  n nodes,  the 


inforaation  for  the  next  level  (^)  nodes  is  produced  and 


kept  in  aeacry  in  order  to  generate  the  inforaation  fcr  the 

next  level,  and  so  on.  Another  approach  is  to  search  the 
graph  depth  first,  keeping  in  storage  cne  only  node  froa 
each  level.  Hhen  the  generated  node  has  its  last  digit 
equal  tc  n , the  procedure  backtracks  in  the  previous  level 
and  the  graph  is  searched  again  all  the  way  down  until  a 
node  with  a last  digit  of  n is  generated.  As  an  exaaple,  in 
a problea  with  5 variables  the  nodes  0,1,12,123,1234  and 
12345  are  first  generated,  and  then  the  procedure  backtracks 
and  replaces  in  level  4 node  1234  with  the  next  successor  of 
node  123,  node  1235.  Since  again  the  last  digit  is  equal  to 
n,  the  procedure  goes  back  two  levels  and  froa  node  12 
generates  nodes  124  and  1245;  froa  here  two  levels  back 
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' I again  and  frca  node  12  generates  now  node  125^  and  so  on. 


In  the  first  approach*  the  nuaber  of  nodes  which  aust  be 
kept  in  storage  changes  froa  level  to  level  and  its  aaziaaa 
n 

value  is  < ) where  l*n/2  for  n even  or  1=  n/2±1/2  for  n 

odd.  In  the  second  approach  the  nuaber  is  constant  and  it 


is  equal  to  n-t-l . lable  4 gives  an  indication  of  these 
nuabers. 


n 

5 

6 

7 

10 

15 

20 

n 

' 1 » 

10 

20 

35 

252 

6435 

184756 

n+1 

6 

7 

8 

11 

16 

21 

. > 

TABLE  4 

■S 


In  the  FORTRAN  prograa  of  the  appendix,  the  second 
I approach  has  been  used. 


« ' « 


r 


I 


The  coaputeE  pcogcaa  for  the  iapleaetatioa  of  the 
algcritha  has  been  sritten  in  FOBIAM,  the  aost  popular 
language  for  the  Operation  Researchers.  The  algoritba  for 
this  prograa  is  nearly  that  described  in  the  previous  [ 

section;  the  sets  are  iapleaented  as  arrays  which  are 
searched  and  updated  when  required.  & block  structured  | 

language  such  as  PASCAL  or  PLI,  peraiting  the  use  of  sets 
and  array  coaparisons,  would  allow  a clearer,  more  faithful  j 

representation.  | 

} 

The  prograa  consists  of  the  main  prcgraa  anc  the  | 

i 

subroutine  CHILD.  In  the  aain  prograa,  the  solution  0 is  ) 

first  tested  for  feasibility.  If  it  is  feasible,  the  next 

I 

node  is  generated  froa  graph  A (node  1)  , and  the  value  of  | 

the  solution  0 is  added  to  the  set  of  efficient  solutions  [ 

(SES) . If  the  solution  0 is  not  feasible,  its  first  j 

successor  node  is  generated  froa  graph  B (i.e.  the  node  j 

which  represents  the  digit  ^*^) . This  is  the  initial  step. 

Froa  now  on  the  recursive  equations  (4)  can  be  used  since  | 

the  values  of  Zo  and  Zc  in  level  0 are  both  zero.  After  the  t 

successor  to  node  0 has  been  generated  either  in  graph  A or  ^ 

in  graph  B,  the  subroutine  CHILD  is  called  tc  test  this  node  | 

and  to  return  to  the  aain  prograa  the  order  to  generate 
(I6MBT>1)  or  to  not  generate  (I6NBT-0)  its  successors  I 

(children) . Defending  on  the  value  of  the  paraaeter  IGNRT,  | 

the  aain  prograa  searches,  always  depth  first,  the  graphs  j 

and  generates  the  next  solution  to  be  tested  by  the  I 

subroutine  CHILD.  Froa  the  aain  prograa  twc  paraaeters  are  ! 
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passing  to  the  subroutine.  The  paraseter  j corresponds  to 
the  level  of  the  node  fron  which  the  exasined  node  has  been 
generated,  and  the  parameter  m indicates  the  last  digit  in 
the  examined  node. 

The  subroutine  CHILD  works  as  follows.  First  the  values 
cf  the  objective  function  for  the  examined  node  are 
calculated.  If  the  SES  is  empty,  the  feasibility  of  the  node 
is  examined  and  if  it  is  feasible  IGNBT~0  is  returned. 
Otherwise  the  order  to  generate  the  child  is  given  tc  the 
main  program.  The  subroutine  continues  to  test  for 
feasibility,  until  finding  the  first  feasible  solution. 
This  solution  is  added  to  SES  and  frcm  that  point  the 
program  tests  first  if  the  examined  node  is  dominated  by 
some  node  in  SES  and  then,  if  it  is  net,  it  tests  its 
feasibility.  If  a solution  is  not  dominated  and  it  is 
feasible,  then  it  is  added  in  to  SES  and  the  solutions  which 
are  dominated  by  the  new  soluticn  are  removed  from  SES.  The 
variable  MSFS  keeps  track  cf  the  number  of  solutions  which 
are  in  SES.  The  values  of  the  variables  for  the  solutions 
in  SES  are  stored  in  the  array  XI,  so  that  the  program  is 
able  to  print  both  the  values  of  the  objective  functions  and 
the  solutions  to  which  they  correspond. 

Using  this  program  the  example  problem  has  been  solved 
and  the  nodes  which  have  been  tested  are  printed,  sc  that 
the  reader  tc  can  compare  these  results  with  the  ones 
obtained  from  Table  3.  The  only  difference  here  is  that 
node  45  has  been  tested;  in  Table  3 the  fact  that  node  5 was 
bounded,  eliminated  the  need  to  test  the  node  45.  The 
reason  is  that  as  long  as  the  graphs  are  tested  depth  first, 
after  node  4 has  been  tested  node  45  is  generated  first  and 
then  node  5. 

t 

It  must  be  noticed  that  this  program  can  be  used  to 
solve  one  objective  function  problems  also;  the  variable  NO, 
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1 

¥ 

i 


i 

i 

j 


I 

I 


which  corresponds  to  the  nuober  of  the  objective  functions, 
is  given  the  value  1. 

The  results  from  17  problems  run  on  the  lEH  360  computer 
of  h.P.S.  have  been  summarized  in  Table  5. 


1 

i number 

1 number  of  | 

cpu 

1 

n 

1 

na 

I no 

1 nc 

1 of  sins 

1 efficient  | 

tine 

1 

i 

1 tested 

1 solutions  1 

sec. 

.1 

5 

1 

2 

1 2 

1 3 

1 14 

1 2 1 

0.24 

1 

6 

1 

0 

1 3 

i 4 

1 21 

1 1 I 

0.33  I 

1 

1 

2 

i 2 

i 4 

i 64 

I 0 I 

0.  66 

« j 

1 

2 

1 3 

1 4 

1 64 

I 0 I 

0.68 

1 

1 

2 

1 4 

i 4 

1 58 

I 5 I 

0.80 

J 

J 

3 

i 3 

1 6 

1 58 

i 3 I 

0.64 

1 

9 

1 

0 

i 3 

i 4 

I 163 

I 2 I 

1.  SI 

1 

1 

3 

i 3 

1 5 

1 346 

I 4 I 

3.67 

1 

1 

3 

1 2 

1 4 

J 377 

( 5 i 

3.39 

1 

3 

1 3 

1 4 

1 97 

i 4 I 

1.  14 

1 

1 

3 

1 

1 4 

1 2 53 

i 13  1 

3.68 

1 

10 

J 

0 

1 1 

1 7 

1 388 

1 1 I 

3.43 

1 

1 

0 

1 3 

1 4 

1 175 

< 3 I 

2.42 

J 

1 

3 

1 2 

I 4 

1 224 

I 2 I 

2.  23 

1 

1 

3 

1 3 

1 4 

i 2 32 

I 8 i 

2.  74 

i 

1 

4 

1 3 

1 4 

1 331 

I 2 I 

4.  19 

1 

25 

1 

5 

i 2 

1 3 

1 747 

I 7 I 

2.53 

1 

1 

5 

1 2 

1 3 

1 12049 

i 7 I 

14.72 

1 

' 

TABLE  5 
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THIS  PROGRAM  Sa  VES  A MULTI  OR  ONE-CBJECTIVE  FUNCTION 
INTEGER  ZERO-CNE  MINIMIZATION  PROELEMt  FORMALIZED  IN 
THE  CANONICAL  FORM. 


CREATOR:  AGGELOS  C.  SIMOPOULOS 
MAJOR  HELLENIC  ARMY 
DECEMBER  1977 


THE  FOLLOWING  NOTATION  HAS  BEEN  USEC: 

N s NUMBER  OF  VARIABLES  IN  THE  PRCELEM 
NA  = NUMBER  OF  VARIABLES  ASSOCIATED  WITH  AT  LEAST  ONE 
NEGATIVE  COEFFICIENT 
NB  * N“NA 

NO  * NUMBER  OF  OBJECTIVE  FUNCTIONS 
NC  » NUMBER  OF  CONSTRAINTS 

ZO(N+1,NO)  » ARRAY  CONTAINING  THE  VALUES  CF  THE 

CBJ.  FUNCTIONS  IN  EACH  LEVEL 
ZCtN+l,NC}  = ARRAY  CONTAINING  THE  VALUES  OF  THE 

CONSTRAINTS  IN  EACH  LEVEL 
SES(**iNO)  = SET  OF  AT  MOST  **  EFFICIENT  SLNS 
IFLSES  » FLAG;  INDICATES  IF  SES  IS  EMPTY  OR  NCT 

IFLAGl  = FLAG;  INDICATES  IF  SLN  C IS  FEASIBLE  OR  NOT 

NSES  = CURRENT  NUMBER  CF  SLNS  IN  SES 

>!(♦♦, N)  = CONTAINS  THE  VALUES  CF  THE  VARIABLES  FOR 
♦♦  SOLUTIONS  IN  SES 
NXl  = CURRENT  NUMBER  OF  SLNS  IN  XI 

LCA(NA+li  = LAST  DIGIT  OF  A NODE  IN  GRAPH  "A"  FOR 
EACH  LEVEL  CF  THE  GRAPH. 

LCB(NB+1)  » AS  ABOVE  FOR  GRAPH  "B" 

MA  * CURRENT  LAST  DIGIT  IN  A NCDE 
LA  a CURRENT  LEVEL  CF  GRAPH  "A” 

LB  a CURRENT  LEVEL  CF  GRAPH  "B" 


C(NC,NJ  a MATRIX  OF  COEFFICIENTS  FOR  THE  0.  FUNCTIONS 
AINCtNI  a MATRIX  OF  COEFFICIENTS  FOR  THE  CONSTRAINTS 
E(NC)  a the  B vector  FOR  THE  CONSTRAINT  ECUATIQNS 


THE  FOLLOWING  PROGRAM  HAS  BEEN  ARRANGED  TO  SOLVE 
PROBLEMS  UP  TC  25  VARIABLESt  5 OBJECTIVE  FUNCTIONS 
AND  10  CONSTRAINT  EQUATIONS.  DEPENDING  ON  THE 
NUMBER  OF  THE  VARIABLES  N,  FORMATS  1201  AND  1202 
MUST  BE  CHANGED  EACH  TIME  TO  NF5.0  AND  (N+1)F5.0 
CORRESPONDINGLY  TO  GIVE  SUITABLE  OUTPUT. 

IF  AN  ERROR  OCCURS  AND  THE  USER  HAS  CHECKED  TO  BE 
SURE  THAT  THE  INPUT  DATA  IS  ACCORDING  TC  THE  FORMAT 
OF  THE  PROGRAM,  THEN  THE  FIRST  DIMENSION  CF  THE 
ARRAYS  XI  AND  SES  MUST  BE  INCREASED. 

FOR  BIGGER  PROBLEMS  THE  DECLARATIONS  AND  FORMATS 
MUST  BE  CHANGED  ACCORDING  TO  THE  NOTATION  GIVEN 
ABCVE. 
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0LCCK  DATA 
INTEGER  X(25) 

DIMENSION  ZO(26t5)f ZC(26f 10) 

COMMON  /S/  XfIFLSESf  IFLAGlfNXltNSES»MA,LA,LE,ZO,ZC 
DATA  X.IFLSESt IFLAGl »NX1 , NSES»MAt LA* LE » ZCt ZC 

1 /29*0»3*1»390#0./ 

END 


INTEGER 

DIMENSION 

L 

COMMON  /S/ 
COMMON  /SI/ 


X(  25)  tXK  150«25),LOA(6)  tLDB(26) 

Z0(26*5) *Z0(26,  10) *0  (5*25) *A( 10«25)*B( 10) t 
SES(  150,5) 

X,IFLSES, IFLAGl ,NX1 , NSES ,MA , LA , L3 , ZO, ZO 
SES,N*N0*NC,C,A*B*X1 


READ  (5,1000  N,NA,NO,NO 
READ  (5,1001)  ( (C(I,J),J-1,N) ,I»1,N0) 
READ  (5,1001)  ( (A(I, J),J»1,N),I»1,NC) 
READ  (5,1002)  ( B( I ) , I *1 ,NC) 


WRITE(6,110C) 

WRiTE(6,1201) 

taRITE(6,1101) 

WRITE(6,1202) 

WFITE(6,il05) 


( (C(I,J),J*1,N),I=1,N0) 
((A(I, J),J=1,N},B(I) ,I»1,NC) 


1000 

1001 

1002 

1201 

1202 

HOC 


1102 

1103 

1104 

1105 

1106 
C 

C 

C 

C 

C 

C 

C 

c 

c 


FCRMAT  (/6X, 25(13)) 

FCRMAT  (415) 

FCRMAT  (10F8.3) 

FCRMAT  (10F8.3) 

FCRMAT  (5F5.0) 

FCRMAT  (6F5.0) 

FCRMAT  (//i2X,' COEFFICIENT  MATRIX  FOR  THE', 
i IX,'  OBJECTIVE  FUNCTIONS' //) 

1101  FCRMAT  (//12X, 'COEFFICIENT  MATRIX  FOR  THE', 

' CONSTRAINTS'//) 

FCRMAT  (////12X, 'THE  SLNS  FOR  THE  ABOVE  PROELEM  ARE') 
FCRMAT  (//5X,'VALUE:' ,5F8.2) 

FCRMAT  (/12X,'S0LLTICN:' ,2514) 

FCRMAT  (///12X,'THE  FCLLOWING  SLNS  HAVE  BEEN  TESTED') 
FCRMAT  (///12X,*N0  FEASIBLE  SOLUTICN  WAS  FOLND* ) 


MISCELLANEOUS 


(NA.EQ.O)  NA«1 
(NA.EQ.N)  NA»N-1 
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^ ^ i 


E ■ 


t . 

I * 

r ■ 


. i 


I 


c 

c 

c 


2 

C 

C 

c 

c 

c 


c 

c 

IS 

c 

c 

c 

c 

c 

20 


22 


c 

c 

c 

c 

c 

120 


130 

135 


138 

C 

bs 


TEST  SOLUTION  3 FIRST 


WRITE  (6»12C0)  (X(I)tI»liN} 

00  2 I«1.NC 

IF  (e(I).GT.O.)  GO  TO  20 
CONTINUE 


OTHERWISE  SOLUTION  0 IS  FEASIBLE 


00  3 I^ltN 
Xl(ltl)»0 
CONTINUE 
N>1>1 

00  4 I^ltNO 
SES(l»n<0. 
CONTINUE 
NSES*1 
IFLSES’l 


IFLAG1=1 

LAal 


GENERATE  NOOES  FROM  GRAPH  A 


CC  30  JA-LAyNA 
P0=NA+1 

IF  ( IFLAGl  .NE.l ) GO  TO  120 
L8*JA+1 

00  22  II-MBtN 
X(ir J=0 
CONTINUE 
X(MA)=1 

WRITE(6»1200)  (X(I}.I»lyN) 
CALL  CHILO  ( JAyMAyl 6NRT) 
IF  ( IGNRT.EC.O)  GO  TO  25 


GENERATE  NOCES  FRCH  GRAPH  8 


00  130  JB=LBtN 
X ( MB ) * 1 

WRITE?6yl200)  (X(n,I»l,N) 
CALL  CHILO( JB.MfiyIGNRT) 

IF  (IGNRT.EQ.O)  GC  TO  13S 
L0B( je)-MB 


IF  (MB.EQ.N) 
M8=MB+1 
CCNTISUE 
LB*  J 0-1 

IF  ( IFLAGl. EQ.O)  L8*LB>1 
IF  { LB.EC.JA)  GC  TO  25 
IF  ( IFLAGl. EQ.O)  LB*LB-1 
LLB*LOB(LB) 

00  138  K*LL6yN 
X(K  )*0 
CONTINUE 
MB*LLB+1 
GO  TC  120 


GO  TC  135 


IF  (MB.EQ.N) 
X(MB)*0 
LB*JB 
MB-MB+1 
GO  TO  12C 


GO  TO  135 
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f 

i 


c 

c 

c 

c 

25 


30 

35 


C 

c 

c 

c 

c 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 


c 

c 

c 

c 

c 


RETURN  TQ  GRAPH  A 


IF  (IFLAGl  .NE.l)  GO  TO  19 
ICA( JA)«HA 

IF  (MA.EQ.  NA)  GO  TO  35 

PA-HA-fl 

CONTINUE 

LA»JA-1 

IF  (LA.EG.O)  GO  TO  80 
LLA«LOA(LA) 

CO  36  KA-LLA«NA 
X(KA)=0 
CONTINUE 
HA»LLAn 
C-C  TO  20 


PRINT  OUT  THE  SOLUTIONS 


80 

IF  (NSES.EQ.O) 

WFITE(6,1102) 

DO  85  1*1, 

WRITE(6 

WRITE(6 

85 

CONTINUE 

GC  TO  98 

95 

UFITE(6,1106) 

98 

STOP 

ENC 

GO  TO  95 


(X1(I» IJ)tIJ«ltN) 


THE  SUBROUTINE  CHILD  TESTS  A SOLUTICN  FOR  DOMINATION 
AND  FEASIBILITY  AND  RETURNS  TO  THE  MAIN  PROGRAM  THE 
CRCER  TO  GENERATE  OR  NCT  GENERATE  THE  SUCCESSORS  OF 
THE  CQRRESPQNCING  NODE 


SUBROUTINE  CHILD  (J.MtlGNRT) 


-MILU  » I 

INTEGER  X(25l»Xl(i50.25) 

DIMENSION  ZO(26t5) yZC(26t 10) «C(5»25) ,A(10»25)tB< 10)t 
1 SES(l50i5) 

CCMMON  /S/  XtlFLSESt  IFLAGlfNXltNSEStMAfLA,LB»ZQTZC 
COMMON  /SI/  SEStN,N0tNCtCfA,6TXi 


202 

C 


DO  202  I>lyNO 

ZO(jnfI)«ZO(JyI)K(I»M) 

CONTINUE 


IF  (IFLSES.EQ.O)  GO  TO  210 
TEST  IF  THE  SLN  IS  DCMINATEO 


:0  5 K-ltr 
ICNTR*6 


DO  204  1*1, NO 


204 


205 

C 


IF  {Z0(J+1,I).LT.SES(K,I))  GO  TC  205 
IF  (Z0( J+1,I} .EQ.SES(K,I) ) ICNTP* ICNTR+1 
CONTINUE 

IF  (ICNTR.EQ.NO)  GO  TO  205 
GO  TO  220 
CONTINUE 
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c 

C CHECK  FOR  FEASIBILITY 
C 
C 

210 


212 


00  212  I-ltNC 

ZC(  J<«>ltI}«ZC(JfI)-«>A(IrM) 
CONTINUE 


215 

C 


00  JjJjNC 


c 

c 

c 

c 

c 


CONTINUE 
IF  (IFLSES.EQ.0)  GC  TC  315 


J-t>l»I).LT.B(  in  GC  TO  230 


UPDATE  SES  AND  XI 


302 

C 


N>1«NSES>1 

DO  302  I«l«N 

Xl(NXl,n*X(I) 

CONTINUE 


C 

C 

C 

C 


ELIMINATE  CQHINATED  SLNS  FRON  SES  AKC  XI 
K»1 


DO  310  I«1.NSES 
■ ITR«0 


305 


ICN'  _ 

DO  305  II-lfNC 
CONTINUE 

IF  (ICNTR.EQ.NO)  GO  TO  306 
GO  TO  310 


C 

C CTHERtilSE  KEEP  THE  SLNS  IN  SES  AND  IN  XI 
C 

306 


308 


DO  308  IK«1,NC 

.iK»«SES( ItIKJ 


305 

310 


312 

313 


SES(K 
CONTINUE 
DO  309  IN>ltN 

XI  (Kf  IN)«X1(  If  IN) 
CONTINUE 

K-K41 

CONTINUE 
DO  312  lK>lfNC 

SES(Kf IK)«Z0( J41,IK) 
CONTINUE 
DO  313  IN-l.N 

Xl(KfIN)«Xl(NXlf IN) 
CONTINUE 
NSES«K 


315 

317 


GC  TO  220 

317  lK«ltNG 
SESdf  IK)>ZO(  J4lf  IK) 


DO 


318 


i 


CONTINUE 
N<ES«l 
IFLSES«1 

DO  318  1>1  ,N 

Xl(  itiUxd) 
CONTINUE 
NX1«1 


220 

30 


IGNRT-0 

RETURN 


IGNRT«1 

RETURN 


0 TC  306 
CNTR»ICNTR4l 


COEFFITIENT  FATRIX  FOR  THE  CEJECTIVE  FUNCTIONS 


COEFFICIENT  MATRIX  FOR  THE  CONSTRAINTS 


1. 

1. 

1. 


k 


1: 


THE  FOLLOWING  SLNS  HAVE  BEEN  TESTED 
C 0 0 0 0 

C 0 1 0 0 

C 0 0 1 0 

C 0 0 1 1 

C 0 0 0 1 

1 0 0 0 0 

10  10  0 
1 0 C 1 0 

1 0 0 0 1 

1 1 C 0 0 

C 1 0 0 0 

C 1 1 0 0 

C 1 0 1 0 

C 1 C 0 1 


THE  SLNS  FOR  THE  ABOVE  PROBLEM  ARE 

VALLE:  2.00  2.C0 

solution:  C 0 1 0 C 

VALLE:  3.0C  -3.00 

solution:  1 1 0 0 C 
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